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Abstract. We give a brief review of violations of the fluctuation-dissipation theorem 
(FDT) in out-of-equilibrium systems; in mean field scenarios the corresponding 
fluctuation-dissipation (FD) plots can, in the limit of long times, be used to define 
a effective temperature T^s that shares many properties of the true thermodynamic 
temperature T. We discuss carefully how correlation and response functions need to 
be represented to obtain meaningful limiting FD plots in non-mean field systems. A 
minimum requirement on the resulting effective temperatures is that they should be 
independent of the observable whose correlator and response are being considered; 
we show for two simple models with glassy dynamics (Bouchaud's trap model and 
the Glauber-Ising chain at zero temperature) that this is generically not the case. 
Consequences for the wider applicability of effective temperatures derived from FD 
relations are discussed; one intriguing possibility is that at least the limit of the FDT 
violation factor for well separated times may generically be observable-independent 
and so could yield a meaningful Teff. 



PACS numbers: 05.20.-y; 05.40.-a; 05.70.Ln; 64.70.Pf 

1. Introduction: Fluctuation-dissipation relations and effective 
temperatures 

One of the core ideas of statistical mechanics is that equilibrium states can be accurately 
described in terms of only a small number of thermodynamic variables, such as 
temperature and pressure. For glassy systems, which can remain far from equilibrium 
on very long time scales, no similar simplification exists a priori; the whole past history 
of a sample is in principle required to specify its state at a given time. This complexity 
makes the theoretical analysis of glasses very awkward, and one is driven instead to look 
for a description of out-of-equilibrium states in terms of a few effective thermodynamic 
parameters. The focus of the present paper is on one such parameter, the effective 

§ Email: peter.sollich@kcl.ac.uk 

II Present address: Department of Physics and Astronomy & Polymer IRC, University of Leeds, Leeds, 
LS2 9JT, U.K. 

^ Email: peter.mayer@kcl.ac.uk 



Effective temperatures in non-mean field systems 



2 



temperature; this can be defined on the basis of fluctuation-dissipation (FD) relations 
between correlation and response functions and has proved to be very fruitful in mean 
field systems 

The use of FD relations to quantify the out-of-equilibrium dynamics in glassy 
systems is motivated by the occurrence of ageing 0: The time scale of response to an 
external perturbation increases with the age (time since preparation) t^ of the system. 
As a consequence, time translational invariance (TTI) and the equilibrium fluctuation- 
dissipation theorem (FDT) relating correlation and response functions break down. 
To quantify this, consider the autocorrelation function for a generic observable m of a 
system, defined as 

tw) = (m(t)m(t^)) - (m(t)) (m(t„)) (1) 
The associated 'impulse response' function can be defined as 

6{m{t)) 



R(t, t^ 



6h{t^ 



h=0 

and gives the linear response of m{t) to a small impulse in its conjugate field h at time 
tw An equivalent way of characterizing the linear response is via the 'step response' 
function 

X(t,tw)= fdt'R{t,t') (2) 

J t-w 

which tells us how m responds to a small step h(t) = hQ{t — t^) in the field. In spin 
glasses, this response function would be called the zero-field cooled response because 
the field is only switched on a certain 'waiting time' t„ after preparation of the system 
{e.g. by cooling) at time t = 0. 

Now, in equilibrium, C(t,t^) = C(t — t„) by TTI (similarly for R and x), and the 
FDT reads 

d \ d 

-X{t-t„) = R{t,t^) = -—C{t-t^) (3) 



dt^-^' " ' Tdt,, 

with T the thermodynamic temperature (we set = 1). A parametric 'FD plot' of x 
vs. C is thus a straight line of slope — 1/T. In the ageing case, the violation of FDT can 
be measured by an FD ratio, X{t,t„), defined through 0, |^ 

^ -X(t, tw) = R{t, tw) = ^%^T^C(t, t^) (4) 



dt^-^' ' ^ ' T dt,. 

In equilibrium, due to TTI, the derivatives d/dt^ in the FDT (|^) could equally be 
replaced by —d/dt; one could therefore argue that a similar replacement could be made 
in (^), leading to a different definition of X{t,t^). However, the latter would make 
rather less sense, since in a situation without TTI only the t^-derivative of x(t,t^) is 
directly related to the impulse response R{t,t^). Adopting therefore the definition (^, 
one sees that values of X different from unity mark a violation of FDT. In glasses, 
these can persist even in the limit of long times, indicating strongly non-equilibrium 
behavior even though one-time observables of the system — such as entropy and average 
energy — may have already settled to essentially stationary values. 
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Remarkably, the FD ratio for several mean field models 0, ^] assumes a special form 
at long times: Taking t^ ^ oo and t ^ oo at constant C = C{t,t^), X{t,t^) X{C) 
becomes a (nontrivial) function of the single argument^ C. If the equal-time correlator 
C{t,t) also approaches a constant Cq for t oo, it follows that 

X(t,tw) = ;^r° dCX{C) (5) 

Graphically, this limiting non-equilibrium FD relation is obtained by plotting x vs C 
for increasingly large times; from the slope —X{C)/T of the limit plot, an effective 
temperature ^ can be defined as Teff(C) = T/X{C). 

In the most general ageing scenario, a system displays dynamics on several 
characteristic time scales, one of which may remain finite as tw ^ cxd, while the others 
diverge with t^- If, due to their different functional dependence on t^, these time scales 
become infinitely separated as tw — ^ cc, they form a set of distinct 'time sectors'; in 
mean field. Teg (C) can then be shown to be constant within each such sector p. In the 
short time sector {t—t^ = 0(1)), where C(t,t^) decays from Cq to some plateau value, 
one generically has quasi-equilibrium with Tes = T, giving an initial straight line with 
slope — 1/T in the FD plot. The further decay of C (on ageing time scales t—t„ that grow 
with t^) gives rise to one of three characteristic shapes: (i) In models which statically 
show one step replica symmetry breaking (RSB), e.g. the spherical p-spin model 0], 
there is only one ageing time sector and the FD plot exhibits a second straight line, 
with Tcs > T. (ii) In models of coarsening and domain growth, e.g. the 0{n) model at 
large n, this second straight line is flat, and hence Tgfj = oo |]^. (iii) In models with 
an infinite hierarchy of time sectors (and infinite step RSB in the statics, e.g. the SK 
model) the FD plot is instead a continuous curve p. 

Tgfj has been interpreted as a time scale dependent non-equilibrium temperature, 
and within mean field has been shown to display many of the properties associated with 
a thermodynamic temperature For example (within a given time sector), it is the 
reading which would be shown by a thermometer tuned to respond on that time scale. 
Furthermore — and of crucial importance to its interpretation as a temperature — it is 
independent of the observable^ m used to construct the FD plot 

While the above picture is well established in mean field, its status in non-mean 
field models is less obvious. To check its validity, one must demonstrate a) that a 

This definition is not quite as stringent as requiring that X{t,t^) be a function of C{t,t^) only, 
at long times iw and t. For some coarsening models at criticality one finds, for example, that 
X{t,t^) — f{t'^C{t,ty^)) for long times, in terms of a scaling function /(•) and a positive exponent 
a Q; X is therefore not a function of C alone. Nevertheless, making t and tw large at constant C, 
one finds a nontrivial function X{C) which in this case is discontinuous, X{C) = /(oo) for C > and 
X{C) = /(O) for C 0. 

* This applies to a wide class of observables, but has not to our knowledge been established in complete 
generality. The restriction arises because the observables that are commonly used in mean field models 
are defined in terms of random fields or couplings, and are therefore uncorrelated with the energy of the 
system. It is possible that for other variables, which are correlated with the energy, different effective 
temperatures would be obtained (S Franz and F Ritort, private communication). This effect is observed 
explicitly in oscillator models with Monte Carlo dynamics . 
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limiting FD plot exists, b) that it gives effective temperatures that are independent of 
the observable-field pair used to calculate C and and c) that the effective temperature 
is constant within a given time sector, i.e. the same for all time scales differing only 
by factors of order unity. One then expects FD plots similar to those for mean field 
systems, composed either of a set of straight lines, or of a continuous curve in the case 
of an infinite number of distinct ageing time sectors. Of course, full independence of the 
effective temperature from the observable considered may be too strong a requirement 
for non-mean field systems, but one would at least expect independence within a large 
class of sufficiently 'neutral' observables; we return to this point in Sec. 

Encouragingly, molecular dynamics (MD) and Monte Carlo (MC) simulations of 
binary Lennard- Jones mixtures [^, as well as MC simulations of frustrated lattice 
gases [|lOl (which loosely model structural glasses, whose phenomenology is similar to 
that of the p-spin model) show a limiting plot of type (i). Oscillator models with MC 
dynamics show similar plots, but the initial short-time regime where one has effective 



equilibrium is absent and one finds only a single straight line [11 1. MC simulations of 



the Ising model with conserved and non-conserved order parameter in dimension d = 2 
and 3 show a plot of type (ii) [1^. And MC simulations of the Edwards- Anderson 



model in d = 3 and 4 |]T3[ give a plot of type (iii). The majority of existing studies, 
however, do not show that Teg is independent of observable (notable exceptions being 
Refs. |10]) since they consider just one observable-field pair. In the Backgammon 
model (see e.g. Ref. Jl^ for a summary of results), several observables have recently 
been considered and differences in the FD ratios found; however, these differences only 
occur in sub-leading terms which decay to zero at long times, and so it is not clear how 
seriously this observable-dependence needs to be taken. 

Below, we therefore give some results for two models for which FD relations can be 
calculated for a broad range of observables; our aim is to clarify the status of FD-derived 
effective temperatures in non-mean field systems. First, however, we pause briefly to 
discuss the appropriate representation of FD plots in such systems. As explained above, 
for mean field systems the existence of a limiting relation (|^) between response x 
correlation C ensures that parametric plots of x versus C converge, for long times, to 
a limiting FD plot whose negative slope directly gives X{C)/T. Eq. (^ implies that 
the plots can be produced either with t as the curve parameter, holding the earlier 
time fixed, or vice versa. Since it is normally easier to fix the time at which a 
field is switched on and measure the resulting response as a function of t — rather than 
observing the response at fixed t to fields switched on at a range of earlier times — the 
first version is the one that is conventionally used 0, |^. 

In general, however, it is clear that the definition (|) ensures a slope of —X{t, t^)/T 
for a parametric x-C plot only if is used as the parameter, with t being fixed. This 
issue becomes particularly important when — for the particular observable chosen — 
the equal-time correlator C{t,t) diverges or decreases to zero for t — > oo, instead of 
approaching a constant value. A series of 'raw' FD plots at increasing t will then either 
grow or shrink in scale indefinitely, preventing the existence of a limiting FD plot for 
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t ^ oo. It therefore appears sensible to normalize the FD plot. The normalization 
factor must be the same for the x-axis and the C-axis and must be independent of t^, 
in order to preserve the connection between X and the slope of the plot. To 'attach' 
the plot to a definite point on the C-axis, it is thus sensible to normalize by C{t,t), 
plotting x{t,t^) = x{t,tw)/C{t,t) vs C{t,tw) = C{t,tvj)/C{t,t). In this representation, 
a limiting FD plot may be approached as t is increased. If this is the case, then X(t, t„) 
converges to a nontrivial limit for large times, but now at a constant value of the 
normalized correlator C (rather than of C itself). Also, if a limiting plot exists, it can 
be obtained from plots of x vs C that use either t or t^ as the curve parameter; the 
normalization still has to be performed with C{t,t) and not with C{t^,t^), however, in 
order to retain the link between X and the slope of the limit plot. 

One final benefit of the method of producing FD plots discussed above — normalize 
by C(t,t), and use tw as the curve parameter — is that the plots for the two types of 
response function most frequently considered are trivially related. We have worked above 
with the 'switch-on' (or zero-field cooled) response function @] a common alternative 
is the 'switch-off' (or thermo-remanent) response, defined for a field switched on at time 
and off again at time t„: 

X_ (t, tw) = / dt' R{t, t') = x{t, 0) - x{t, tw) (6) 
Jo 

With tvi being used as the curve parameter, the FD plots of x vs C and X- vs C are 
then trivially related (by a reflection about the horizontal axis and a vertical shift); the 
same remains true after normalization of the curves by C{t,t). If, on the other hand, 
one were to follow the naive approach of using t as the curve parameter — keeping t„ 
fixed — and normalizing by C{t^,t^), then the plots for x and for X- are not related in 
any obvious way and can in fact have very different shapes |^ . 



We end this section with a brief comment on a recent proposal [T^ to construct 
FDT plots from the response and the disconnected correlation function 

Cdis(t,tw) = {m(t)m{t„)) 

as opposed to the connected one (|l|) used above. This can avoid the normalization 
problem in some cases: If, for example, one is dealing with a systems of spins Si = ±1, 
and m is a staggered magnetization with a random field, then Cdis reduces to a spin-spin 
autocorrelation function and so is automatically normalized at equal times, C^isit, t) = 1. 
While apparently convenient, use of Cdis in FD plots appears to us to have little 
theoretical justification. A serious drawback is that FD plots constructed from Cdis are 
no longer invariant under simple affine transformations of the observable, m' = am + b. 
This seems to us to violate the physical intuition that such trivially related observables 
should not be able to reveal significantly different glassy physics. We would therefore 
advocate use of the connected correlation function (|I|) throughout; this ensures that 
affine transformations of observables, as defined above, leave all FD plots invariant (up 
to a trivial overall scale factor of a^). Of course, if the observable is such that it has no 
drift, i.e. (7Ti(t)) = constant, then disconnected or connected correlation functions can 
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be used interchangeably, since they differ only by a trivial additive constant. This may 
be the case for sufficiently 'neutral' observables (or could in fact be part of the criterion 
for deciding whether an observable is neutral); see the discussion in Sec. ^ 

2. Model system 1: Bouchaud's trap model 

Trap models |T^, [TB|, |TP| are obvious candidates for a study of FD violation in systems 
with glassy dynamics. Popular as alternatives to the microscopic spin models discussed 
above, they capture ageing within a simplified single particle description. The simplest 
such model [l^ comprises an ensemble of uncoupled particles exploring a spatially 
unstructured landscape of (free) energy traps by thermal activation. The tops of the 
traps are at a common energy level and their depths E have a 'prior' distribution p{E) 
{E > 0). A particle in a trap of depth E escapes on a time scale t{E) = TQexp{E/T) 
and hops into another trap, the depth of which is drawn at random from p{E). The 
probability, P{E,t), of finding a randomly chosen particle in a trap of depth E at time 
t thus obeys 

id/dt)PiE,t) = -T-\E)P{E,t) + Y{t)p{E) (7) 

in which the first (second) term on the RHS represents hops out of (into) traps of depth 
E, and Y{t) = (r^^(£'))p^^ is the average hopping rate. The solution of (0) for initial 
condition Po{E) is 

P{E, t) = Po(^)e-*/"(^) + p{E) fdt' F(t')e-(*-*')/"(^) (8) 

•J 

from which Y{t) has to be determined self-consistently. For the specific choice of prior 
distribution p{E) ~ exp {—E/Tg), the model shows a glass transition at a temperature 
Tg. This can be seen as follows. At a temperature T, the equilibrium Boltzmann state 
(if it exists) is Peq{E) oc t{E)p{E) oc exp (E/T) exp {—E/Tg). For temperatures T <Tg 
this is unnormalizable, and cannot exist; the lifetime averaged over the prior, (t)^, is 
infinite. Following a quench to T < Tg, the system never reaches a steady state, but 
instead ages: in the limit t„ — oo, P{E, t^) is concentrated entirely in traps of lifetime 
r = 0(tw)- The model thus has just one characteristic time scale, which grows linearly 
with the age t„. In contrast, for T > Tg all relaxation processes occur on time scales 
0{tq). In what follows, we rescale all energies such that Tg = 1, and also set tq = 1. 

To study FD violation one can extend the model as follows |2^: To each trap 
we assign, in addition to its depth E, a value for an (arbitrary) observable m; by 
analogy with spin models we refer to m as magnetization. The trap population is 
then characterized by the joint prior distribution a{m\E)p{E), where a{m\E) is the 
distribution of m across traps of given fixed energy E. We focus on the non-equilibrium 
dynamics after a quench at t = from T = oo to T < 1; the initial condition is thus 
Po{E^m) = a{m\E)p{E). The subsequent evolution is governed by 

{d/dt)P{E, m, t) = -T-\E, m)P{E, m, t) + Y{t)a{m\E)p{E) (9) 
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where the activation times are modified, in the presence of a small field h conjugate to 
m, according to T{E,m) = t{E) exp {mh/T). Other choices of T{E,m) that maintain 
detailed balance are possible [|1^, |21[]; we adopt this particular one because, in the spirit 
of the unperturbed model, it ensures that the jump rate between any two states depends 
only on the initial state, and not the final one. 

Let us outline briefiy the calculation of correlation and response function for the 
observable m. The autocorrelation function is (in the absence of a field, i.e. h = 0) 

C{t,t^) = jdm dm^ dE dE„ mm^ P{E^,m„,t^) x 

X [P{E, m, t\E^, m„, t^) - P{E, m, t)] (10) 

in which the 'propagator' P{E,m,t\E^,m„,tvj) is the probability that a particle with 
magnetization m„ and energy E^^ at time t„ subsequently has m and E at time t. 
The key simplification in the calculation is that this can be related back to the known 
distribution P{E,m,t): 

P{E, m, t\E^, m^, t^) = 6{m - m^)6{E - E^)e~^'-'-^'^^''-^ 

+ tdt'T-\E^)e-^''-'-^/^^''-^P{E,t-t')a{m\E) (11) 

The terms corresponds to the particle not having hopped at all since t^ (first term) 
or having first hopped at t' (second term) into another trap; after a hop the particle 
evolves as if 'reset' to time zero since it selects its new trap from the prior distribution 
ij{m\E)p{E), which describes the initial state of the system. We also have P{E,m,t) = 
a{m\E)P{E,t) since at zero field the unperturbed dynamics is recovered. Substituting 
these relations into (|10D and integrating over m and m^, one obtains an exact integral 
expression for C. This is expressible as the sum of two components: The first depends 
only on the mean rn{E) of the fixed energy distribution (T{m\E) and the second only 
on its variance, A'^{E). This of course makes sense, since the correlation function only 
measures statistics of m up to second order. 



To find the corresponding response function, one can use the relation |20] 



T^xit, tw) = ^C{t, t„) + ^ (m(t)) (m(t.)) (12) 



in which {m{t)) = {m{E))p^^^-^ is the global mean of m. Equation (|T2D is a slight 



generalization of the results of [^, it is exact for any Markov process in which 
the effect of the field on the transition rate between any two states depends only on 
the initial state and not the final one. Substituting the expression for C into ([T2|) , 
and integrating over t„ with the boundary condition x(^w,^w) = 0, one finds a closed 
expression for x- For convenience, we rescale the field h — Th in the following, thus 
absorbing a factor 1/T into the definition of the response function. In this way, the slope 
of the FD plot of x vs C, or of their normalized analogues, becomes —X = —T/T^q 
(= —1 in equilibrium). 

Using the exact results described above, C and x have been calculated for a 
number of different distributions a{m\E)., each specified by given functional forms 
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Figure 1. FD plots oi x vs C for a distribution a{m\E) of variance exp{nE /T) 
(but zero mean) for n = 0.2, 0.0, —0.2; T = 0.3. For each n data are shown for times 
t = 10^, 10^; these are indistinguishable, confirming that the limiting FD plot has been 
attained. Dashed line: The predicted asymptote x^l^C'fori— *(X) and C* — > 1. 
From Ref. [|o). 

of m{E) and A'^{E) p^. Each form effectively corresponds to a distinct physical 
identity of the observable m. For example a{m\E) = 6{m — E) implies m = E, in 
which case h = —ST/T. For simplicity, we confined ourselves to distributions either of 
zero mean (but non-zero variance) or of zero variance (but non-zero mean). Even for 
these simpler observables m, the equal-time correlator C{t,t) can depend strongly on 
t; this is in contrast to the spin models often considered in mean field studies, where 
C{t,t) is automatically normalized. The following plots are therefore of the normalized 
correlation and response functions, C and x', as explained above, these are obtained by 
dividing the raw quantities by the equal-time correlator C{t,t) at the later time. 

For the first class of observables (with zero mean m{E)), we considered a variance 
A'^[E) = e:K.p[En/T) for various values of the exponent n, generalizing the results of 
for n = 0. The equal-time correlator is then C{t, t) = J dE P{E, t) exp (nE/T). For large 
t, its scaling behavior can be deduced from that of P{E,t), and depends on the value 
of n: For n < T — 1, C{t,t) is sensitive only to shallow traps (the population of which 
depletes in time as due to ageing, giving C{t,t) ~ t'^^^) and one correspondingly 
expects the two-time correlator to decay on time scales t — t„ = 0(1), probing only 
quasi-equilibrium behavior. In contrast, for T — 1 < n < T one finds that C{t,t) is 
dominated by the contribution from traps with depths E such that t{E) = 0{t); this 
results in the scaling C{t,t) ~ t". For such n, the two-time correlator will decay on 
ageing time scales t — t^ = 0(t^), and we expect strong violation of equilibrium FDT. 
The numerical results confirm the above predictions. For values of n < T — 1 (not 
shown) the normalized FD plot approaches a straight line of equilibrium slope —1 at 
long times t ^ oo. In the regime T — 1 < n < T (Fig. |I|), equilibrium FDT is strongly 
violated, as expected, but a limiting non- equilibrium FD plot is nevertheless approached 
at long times. This can be proved analytically by showing that C and x share the same 
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Figure 2. FD plots oi x vs C for a distribution with mean exp{nE/2T) (but zero 
variance), for n = 0.2, —0.2; T = 0.3. Curves are shown for times t = 10^, 10^, but are 
indistinguishable except for the zoom-inset on the left hand (upper curve: t = 10^). 
Dashed Une: The predicted asymptote x — ^ ^ C ior t ^ oo, C ^ I. From Ref. [ pO[ . 



scaling variable {t—t^)/t^ in this limit. The slope of each plot varies continuously with 
C. In contrast to mean field, this is not due to an infinite hierarchy of time sectors; 
the variation is in fact continuous across the single time sector t — tw = C'(tw)- More 
seriously, different observables give different plots. One has to be slightly careful here: 
Observable-independence requires that X(t,t^), as a function of the long times t and 
tvj, be independent of the observable considered. The limiting FD plots, however, give 
X as a function of C rather than of t and One thus has to convert the results 
for X derived from different observables back into the time domain before comparing 
them. In our case this means comparing the X{C) values corresponding to a fixed value 
{t — tvj)/t„] such a comparison reveals that the values of X do indeed differ among the 
observables considered. 

In Ref. [^, we also considered observables defined by a fixed energy distribution 
a{m\E) with zero variance and mean rn{E) = exp{En/2T); in this case m is a 
deterministic functions of E. Here too, we found that for n < T — 1 the limiting 
FD plot is of equilibrium form, while for T — 1 < n < T (Fig. ^ a non-equilibrium 
FD plot is approached as t oo, with a shape dependent (here very obviously) on 
the observable m. Note that the response function can now change sign from positive 
to negative for increasing time separation t — t„. This is not untypical for acticated 
dynamics, and can be understood as follows: Suppose for simplicity that n > 0. Then 
rn{E) increases with E and so in equilibrium the application of a positive field conjugate 
to m would increase E and thus m on average, giving the conventional positive response 
which we observe for t close to t„. On the other hand, since the application of the 
field effectively increases the depth of all traps it also slows down the out-of-equilibrium 
relaxation of E and m to larger values, thus giving a negative response for large t — t^. 

Observables with a power law dependence of the mean, m{E) = E^, were also 
considered, with zero variance A'^{E) = as before so that m = m{E) is a deterministic 
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function of E. Then p = 1 corresponds to energy-temperature FD, since the observable 
is energy (m = E) and the conjugate field is, apart from trivial scale factors, the 
temperature. Interestingly, in this case no limiting plot exists. This is because the 
amplitude of the correlator remains finite as t ^ oo, while the response function, for 
any fixed value of C, diverges as Int. 

A final observation to be made from Figs. |I| and ^ is that, whenever a limiting FD 
plot exists for the observables considered here, its slope in the limit of well separated 
times (t ^ tw, corresponding to C{t,t^) — > 0) is always Xqo = 0. 

3. Model system 2: Glauber-Ising chain at zero temperature 

As the second model system, we consider a chain of spins Si with the Ising Hamiltonian 
H = — JI]jSjSj+i. The time evolution is assumed to be given by Glauber dynamics; 
this means that at each time step a spin i is chosen randomly and its state assigned to 
Sj = ±1 with probability proportional to the Boltzmann factor exp{—H/T) for the new 
configuration. Time is scaled so that each spin is updated once, on average, per unit 
time. 

When this system is quenched to a low temperature from a random initial state, 
corresponding to equilibrium at T = oo, the dynamics will progressively line up 
neighbouring spins with each other and alternating domains of up and down spins will 
grow (coarsen) until they reach the equilibrium domain length ^eq ~ exp(2J/T) at a 
time Teq ~ exp(4J/T). In the limit of T — > 0, these scales grow to infinity and the 
systems remains out of equilibrium at any finite time; this will be the limit of interest 
to us. Since the Ising spin chain has a phase transition to an ordered ferromagnetic 
state at T = 0, one can regard the out-of-equilibrium dynamics for T — > as critical 
coarsening Physically, however, this is somewhat different from coarsening at a 
finite-temperature critical point, since in the Ising chain the growing domains are not 
critical but have the structure of either of the ordered phases (spins all up or all down). 

The FD relations for the Glauber-Ising chain have been studied previously for the 
case of the spin-spin autocorrelation (see Ref. |2^; a brief summary of the results can 
be found in Ref. 0). It is a simple matter to extend the calculation to a somewhat 
more general class of observables, defined as 

m = ^eiSi (13) 

i 

where the are quenched random variables with zero mean and translation-invariant 
covariances [eiej] = qj^i; here [. . .] denotes the average over the quenched disorder. The 
corresponding perturbation in the Hamiltonian is —hm, as usual, so that the local field 
acting on each site is hei. With different choices of the correlation function one can then 
interpolate between the two extremes of a uniform field {qk = 1, in which case m is the 
conventional magnetization) and a locally random field (g^ = ^w, the case treated in 
Ref. [0). 
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Because {si{t)) = 0, m as defined by (0) also has an average of zero at all 



times. Using the known results for the two-time spin-spin correlation function (|23 
Cij{t,tv/) = (sj(t)sj(t„)), which can be written as Cij(t,t^) = Cj-i(t,t^) due to 
translation invariance, one then has for the autocorrelation function of m (scaled by 
N to get a quantity of order unity) 

C(t,tw) = ^Zl^i^i9-i(^'^w) = iTTY.^l^l+k] Ck{t,t^) = Y.1kCk{t,t„) (14) 

ij k \-^^ I ) k 

where in the last step we have used that in the limit of an infinitely large system 
the spatial average of e;e/+fc becomes identical to the quenched average [e^e^+fc] = qh- 
The scaled response function of m to its conjugate field is obtained similarly as 
i?(t, tw) = Yl,k Qki^kit, tw), where the translationally invariant response function rj_i(t, t„) 
gives the response of a spin at site i and time t to a local field applied at site j at time t^- 
Using the exact solutions for Ck(t, t„) and rk(t, t„), one can thus evaluate the correlation 
and response functions of m, at least in principle. 

Before giving the results, let us briefly discuss what behaviour one would expect. 
We can define a correlation length of the random fields hi = hei as = YZL-oo Ik- 
This length needs to be compared with the growing spin-spin correlation length ^{t) of 
the coarsening system, which is of the same order as the typical domain length. For 
early times the field is essentially uniform on the lengthscale of the spin correlations 
(^(t) ^ ^h) and so we expect to see behaviour corresponding to the uniform field limit. 
For late times, the field correlations are much more short-ranged than those of the spins 
{C{t) ^ ^/i); across each typical spin domain, the ej change signs many times and we 
expect qualitatively the same effects as for completely uncorrelated fields. The crossover 
between these two regimes takes place when C,(t) ^ C,h- Since the growth of ^{t) is driven 
by the random- walk motion of the domain walls, it scales as ^(t) ~ t^/^, and thus the 
crossover time is of order = ^l- 

To simplify the theoretical analysis, it is convenient to look at the limit where all 
times {t, t„ and Th) are large compared to unity. This then also implies that the field 
correlation length is large, so that we can write qk = q{k/C,h), where q{x) is a smooth 
function of the scaled distance x = k/^h- We leave all details of the calculation to a 
forthcoming publication and display only the results. 

In Fig. |^(a) we show the unnormalized (except for a constant scale factor) FD plot 
for the case of a Lorentzian field correlation function q{x) = 1/(1 + x"^). As before, the 
earher time t„ is used as the curve parameter and t is held constant for each curve. As t 
increases, the amplitude of both correlator and response initially grows, but eventually 
saturates for t^r^. For such long times, the FD plot is identical to that for uncorrelated 
fields 1^ as predicted, with a negative slope X = 1 for t close to tw, corresponding to 
quasi-equilibrium, and X = X^o = 1/2 in the limit of well separated times t ^ t^. For 
shorter times t <^ Th, on the other hand, the FD plot is a straight line of slope X = 1/2 
throughout; this is more easily seen once correlation and response are normalized by 
C(t, t) as in Fig. ||(b). As we will see very shortly, and consistent with our expectations. 
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Figure 3. FD plots for Glauber-Ising chain, for a linear observable (|l^) with fields Cj 
that have a Lorentzian covariance function qk — q{k/£,h) with q{x) ~ 1/(1 + x^). (a) 

1/2 

Raw, unnormalized, response and correlation functions. The constant factor 1/t^ 
only serves to remove the trivial scaling of the amplitude of these functions with the 

1/2 

field correlation length = t^J . Curves are for increasing ratios of t/r/j as shown in 
the legend, (b) FD plots of the normalized correlation and response, C and x- The 
normalization shows that in the limit t/xh <?C 1 one obtains a straight line of negative 
slope X = 1/2, exactly as in the long time limit for uniform fields; see Fig. |[ 



the FD plot in this regime, i.e. before the crossover at t ~ r/j, is identical to that 
obtained for uniform fields. 

Fig. I shows the limiting {t oo) FD plots for a range of field correlation functions 
q{x). At one extreme, we have the case of uniform fields {q{x) = 1), giving a straight 
line of negative slope X = 1/2. At the other extreme, all field correlation functions 
q{x) for which a finite correlation length can be defined have the same limit plot as for 
uncorrelated fields. Because we defined the correlation length above as C,h = J^kL-oo Iky 
the scaled correlation function must obey / dxq{x) < oo for a finite correlation length to 
exist; this requires that q{x) decay faster than x~^ for large x. Correlation functions q{x) 
which decay asymptotically as q{x) ~ a;~" with < n < 1 give intermediate FD plots, 
whose shape is determined only by the asymptotic decay exponent n. Results for several 
values of n between and 1 are shown in Fig. ^; they demonstrate clearly that even 
among the relatively restricted class of observables defined by (|r^) , nontrivially different 
limiting FD plots exist which give different results for the FD ratio X. Intriguingly, 
however, the limit of X for strongly separated times t ^ t^, is Xoo = 1/2 for all 
observables considered here. 

4. Discussion 

Let us summarize our findings. For the trap model, on the basis of exact relations 
for correlation and response functions, we have studied FD relations in the glass phase 
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Figure 4. Limiting FD plots for the Glauber-Ising chain, for different asymptotic 
behaviours q{x) ^ x^" of the field correlation function; the value of n is shown in 
the legend. For n ^ one has the case of uniform fields, which gives a straight line 
of negative slope X — 1/2. The opposite extreme n > 1 corresponds to correlation 
functions with a finite correlation length; these all give the same limit plot, identical 
to that for entirely uncorrelated fields. Values of n between and 1 give nontrivial 
intermediate FD plots. Note that the negative slope X of all plots is, in the limit of 
well separated times, Xoc — 1/2. 



for a broad class of observables m. For observables which indeed probe the ageing 
regime, equihbrium FDT is violated; in most limiting non-equilibrium FD plot 

is approached at long times. In contrast to mean field, this plot depends strongly on 
the observable. It furthermore has a slope which varies continuously across the single 
time sector t—t^ = 0(t„). This shows that in this simple, paradigmatic model of glassy 
dynamics, the mean field concept of an FD-derived effective temperature T^s cannot 
be applied. One may dismiss the trap model as too abstract for this to be of general 
relevance; but by placing traps onto a ci-dimensional lattice and allowing hops only to 
neighbouring traps, a more 'physical' model (of diffusion in the presence of disorder) 
can be obtained. Recent work [jl9[ shows the scaling of correlation functions to be 
unaffected by this modification; since (|1^) would also continue to hold, the results for 
the trap model should be qualitatively unchanged. 

Our results for the Glauber-Ising chain were qualitatively similar. There, we 
restricted ourselves to observables that are linear in the spin variables; again, exact 
results for the correlation and response functions of such observables can be obtained. 
Limiting FD plots always exist, but depend nontrivially on the observable. Except for 
the case where the observable is the magnetization (and the corresponding perturbation 
therefore a uniform field), the slope of the limit plots varies continuously; as in the trap 
model, this variation takes place within a single time sector t— = 0(t^). 

The results above demonstrate that the idea of an effective temperature Tefj derived 



Effective temperatures in non-mean field systems 



14 



from FD relations does not apply straightforwardly to generic non-mean field models. 
Could the concept nevertheless be rescued? One difficulty in the two models considered 
was the non-uniqueness of the limiting FD plots. One could argue that in order to probe 
an inherent Tefj characterizing the non-equilibrium dynamics, the statistical properties of 
the observable must not change significantly across the phase space regions visited during 
ageing. (In coarsening models, similar arguments have been used to exclude observables 



correlated with the order parameter ||12|-) Applying this idea to the trap model, where 
the typical trap depth E increases without bound for t — > oo, a 'neutral' observable 
would presumably require that Tn{E), A'^{E) const, as E oo; with this restriction 
we indeed get a unique limiting FD plot. In the Glauber-Ising chain, the energy E = (H) 
decreases towards the ground state energy Eq during coarsening. A neutral observable 
might then be one whose properties within a microcanonical ensemble at given E do 
not change significantly as E ^ Eq. It is easy to see that, for our observables (0), this 
criterion again produces a limiting FD plot: It only allows field correlation functions 
with a finite correlation length — which all give the same limit plot, namely that for 
un correlated fields — while other observables have a variance which diverges as E ^ Eq. 
The selection of neutral observables may in fact obviate the normalization problem 
for FD plots discussed in detail in this paper: If an observable has a variance C{t,t) 
which either converges to zero or diverges as t increases, then this in itself implies that 
its statistical properties change strongly while the system ages; conversely, for neutral 
observables C{t,t) should approach a finite limit for t — > oo. 

Even if one can, by a judicuous choice of a 'neutral' observable, produce a unique 
limiting FD plot and thus a unique FD ratio X{t,t^), one is still faced with the problem 
that X may — as in our cases — vary continuously across a single time-sector. This almost 
certainly makes X inappropriate for the definition of a meaningful effective temperature 
Tgff = T/X, since two thermometers probing time scales which differ only by factors of 
order unity should be required to measure the same temperature. There is, however, 
the intriguing possibility that the limit of X obtained for long times and large time 
separations, Xoo = limf^^oo linit_»oo X(t, t„) may still correspond to a meaningful T^s- In 
fact, in the two models considered here, X^q was the same for all observables considered, 
as long as limiting FD plots actually existed. These ideas should be explored in future 
work; we are currently in the process of calculating FD plots for second orc?er observables 
such as m = J2i ^i^iSi+i in the Glauber-Ising chain. If the same X^c results, this would 
strongly support the hypothesis that X^o is indeed observable-independent. 

Finally, it is of course likely that there are different classes of non-mean field systems 
as far as the concept of an effective temperature is concerned. One could conjecture 
that FD plots would provide a diagnostic tool here: The appearance of rounded FD 
plots like those shown above may signal that a system does not have a meaningful Teg, 
while limiting FD plots consisting of two straight line segments — as found in simulations 
of some non-mean field models [T^, |T2[ — may indicate that a well-defined effective 
temperature does indeed exist. One of the main tasks ahead must be to develop physical 
criteria that allow us to predict to which of these two classes a given system belongs. It 
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is probable that a strong separation of the dynamics into fast and slow time scales will 



play a role (see e.g. |25|), but this may not be sufficient. The trap model is again 
instructive here: In the glass phase the fast processes (on time scales of order unity) do 
not contribute to the dynamics, and all relevant relaxation processes happen on slow 
time scales of the same order (~ A naive time scale separation argument would then 
suggest that the FD plot should consist of a single straight line, with a non-equilibrium 
slope X 7^ 1, while the exact results presented above show that the actual situation is 
rather more complicated. 
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